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Abstract 

A practical method of calculating weights implementing the Bose - Einstein 
interference effects is presented. The usefulness of the method is demonstrated for 
the events generated by the JETSET/PYTHIA code to describe the UAl proton- 
antiproton data at 630GeV. A good description of data is achieved with a reasonable 
value of the Gaussian width of a two-particle weight factor, which is the only free 
parameter in our calculations. 
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1 Introduction 



Recently various methods of imitating tlie Bose - Einstein interference effects in Monte 
Carlo generators have been discussed [1-3]. This problem became acute with the advent 
of new data on WW production, since quite conflicting estimates of the W mass shifts 
resulting from these effects have been published [4-5]. 

It has been shown that the most commonly used method of momenta shifting (em- 
ployed in the new versions of JETSET [6]) does not reproduce the input shape of the Bose 
- Einstein correlation factor, even after improvements of the original procedure [3]. Thus 
modifying this shape to fit the data we cannot hope to learn much about the space-time 
evolution of the production process. 

On the other hand, the method attributing weights to generated events, much better 
justified theoretically [7], as recently discussed in more detail by Bialas and Krzywicki [8] 
is not easy to implement. The factorial increase of the number of terms to be calculated 
with multiplicity makes impossible the direct application of the method for high energies. 
Only recently Wosiek [9] has indicated a possible scheme shortening such calculations. 
Symmetrizing separately the particles from each hemisphere, as proposed by Haywood 
[1], introduces the unknown bias against the effect for slow particles (in the CM). In 
addition, this procedure does not remove the fundamental difficulties but only shifts them 
to higher multiplicities. Another problem of the weight method is the serious distorsion 
of the multiplicity distribution, since weights enhance the high multiplicity tail [2]. 

Solution to these problems has been proposed recently by Jadach and Zalewski [5], 
who reduce the number of terms replacing the original Bialas - Krzywicki formula by 
an approximation based on the clustering algorithm. The initial average multiplicity is 
restored by rescaling the weights with a simple cV" factor. 

In this note we propose another approximation which seems to be quickly convergent 
to the full formula without prohibitive increase in the computing time. We present the 
method in the next section and then apply it to the description of the UAl data [10]. 
The results are encouraging. We hope to use the same method to clarify the problem of 
the W mass shifts in hadronic decays. 

2 Calculating weights for multiparticle states 

According to the Bialas - Krzywicki prescription, the symmetrization of amplitudes re- 
quired by the Bose - Einstein statistics may be approximated by generating unsym- 
metrized distributions and correcting them a posteriori by the multiplicative weights at- 
tributed to each event. Such a weight is a sum over all permutations of identical particles 
of the products of two particle weight factors WiP(j) calculated for the pair of momenta 
(of i-th particle and the particle, which occupies the i-th place in the permutation P{i})- 

n 

W{n)^ E n^^PW- (1) 
{P{i)} i=i 

Since all factors are positive and wu = 1, the resulting weight is not smaller than one 
(a contribution from identity permutation). One may correct results to keep, e.g., the 
average number of particles fixed; we return to this point later. 
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Since most of the particles detected in experiments are pions, the final weight should 
be actually given by a product of weights calculated separately for positive, negative and 
neutral pions. In fact, the BE interference for neutral particles is not observable (apart 
from the possible effects for hard photons [11]): neutral pions decay before detection, and 
for the resulting photons the effective source size is so big that the BE effects must be 
negligible for momentum differences above a few eV. However, the procedure should not 
change the observable correlations between the numbers of charged and neutral pions. 
Therefore weights for all signs of pions must be taken into account. 

Thus in principle the only arbitrary factor is the function of the difference of two 
momenta Wij{pi — pj). It is natural to try as a first guess the Gaussian function of four- 
momentum difference squared 

Wij^e^P^-P^^'/^"'. (2) 

Of course, different components of momentum difference squared may be multiplied by 
different coefficients, and the shape may be modified. In this note we do not discuss these 
possibilities. Therefore the only parameter is a Gaussian half-width of the distribution a. 

We have tried first to implement this prescription generating by the JETSET/PYTHIA 
program [6,12] the samples of 10^ events of pp collisions at 10 and 30 GeV CM energies. 
For each event for each permutation the product of two-particle weight factors (with 
a =0.14 GeV which corresponds to 1 fm radius of the Fourier transform of the weight 
(2)) is computed and all contributions are added to calculate the weight for the event 
which is used to produce distributions (to be compared with those for all weights equal 
one). 

Unfortunately, for more than ten pions of a given sign the calculations become pro- 
hibitively long. This does not happen at 10 GeV, but already at 30 GeV there are more 
then hundred events with such multiplicities. To get any results we had to exclude them 
from standard weight calculations (attributing to each of them the same value of the 
weight as obtained for the previously generated event). This removes, however, the fiuc- 
tuations which are most interesting for the investigation of very short range correlations. 

Thus here we have separated the sum of all the n! permutations into terms where only 
the permutations which change places of exactly K particles are taken into account: 

^ = (3) 

K 

These terms for X < 6 are 
^(0) ^ 1. ^(1) ^ 0; = Y:7=iEj>Mjy; w^^'^ ^2j:7=iEj>iJ:k>jWijWjkWki; 

^(4) ^ J21-^Y.j>iT,k>jT,l>k['^WijWikWjiWkl + 2WijWiiWjkWkl + 2WikWiiWjkWji + {wuWjkf+ 
{WijWklf + [WikWjif]; 

-ti;(5) = 2 Er=l^ Ej>j Efe>j Ei>ifc Y.m>l[{Wijfwi),WmlWkm + {WikYwjlWmlWjm+ 
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{WiifWjkWjmWkm + {WimfWjkWklWji + {WjkfwuWimWim + {WjlfwikWkmWim+ 

fWikWklWii + {Wkl)'^ + {wirafWijWjkWik + {WkmfWijWjmi + 

WijWjkWklWimWim + WikWjiWkmWjmWil + WilWijWklWjmWkm + WijWikWjiWimWkm+ 
WikWimWjkWjiWim + WuWjiWjkWkmWim + WijWikWklWimWjm + WijWuWimW jkWkm+ 
WijWimWjiWkmWkl + WikWuWjkWimWjm + WikWimWjrn^jlWkl + WuWimWklWjkWjm]- 

(4) 

To check the method we have first compared the results from the full sum of permu- 
tations (1) with the results from the sum (3) cut at = 4 at 10 GeV. Both programs 
give the same results within a few permille accuracy for all investigated distributions. Of 
special interest is the " BE ratio" , defined for the pair of identical pions as a function of 

Jd^Pid^P2p2{pi,P2)S[Q - a/-(pi -P2)2] <n>2 

C2[Q) ~ — , = 7 r . (5) 

fd^Pld^P2Plipi)piip2)S[Q - ^-{Pl-P2f] < - 1) > 

Without weights it is rather fiat and close to one, if we normalize separately the 
numerator and the denominator of Eq. (5) to the same number of entries (which is 
achieved by the second factor in (5)). Including weights produces a maximum at smallest 

with the height about 2 (i.e. one unit above the value at large Q^) and a width 
a' about 0.15. Thus we reproduce satisfactorily the shape assumed for the two-particle 
weight factor. 

At 30 GeV we do not have, as noted above, the results of full symmetrization; for 
174 events of the highest multiphcity the weights of the previously generated events were 
attributed (still, this program requires 10 times more computing time than the program 
with no more than 4 momenta symmetrized!). The multiplicity distribution for two pro- 
grams differs slightly in the tail, although average multiplicities are quite similar: 10% 
and 11% higher than without weights. The peak at low Q'^ in the ratio of distributions 
exceeds slightly 2 and looks similar in both programs. 

After this exercise we started generating events at 630 GeV, the energy of the UAl 
experiment. For the calculations of weights we used the same value of a as before. We did 
not have now the possibility to estimate the results of full symmetrization (the highest 
multiplicities of one sign pions exceeded 40 in the 10^ events sample). Thus we have 
checked first that cutting the series (3) a,t K — 3 and a,t K — A we get quite similar 
shapes of the spectra, although the normalization is significantly different. Including 
the term with K = 5 we change even less all the distributions. Thus we feel that cutting 
the series (3) at K = 5 we get a reliable estimate of the results for spectra from the 
weight method (up to the possible change of normalization). 

The distribution of weights at 630GeV is much broader than at previous energies and 
has a long tail (up to the values of a few hundreds). Consequently, the multiplicity distri- 
bution is significantly changed by the weighting. Since the JETSET/PYTHIA parameters 
were fitted to reproduce inclusive experimental data without weights, the change, e.g.. 
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of the average multiplicity induced by weights should be compensated by the proper re- 
fitting procedures. Instead we have applied a simple method of multiplying weights by 
an extra cV"' factor, where n is the number of pions, and c and V are constants fixed 
by the requirements to restore the original number of events and the original average 
multiphcity. We return later to the details of this procedure. Such a rescaling of weights 
does not change significantly the shape of distributions. The BE ratio refiects mainly 
the assumed shape of the two-particle weight (plus 1): for larger a it is wider and starts 
to increase above 2 for smallest Q^. 

The procedure seems to produce too high a value of the BE ratio for smallest Q'^. As 
already noted, it is about twice the value for small Q^, whereas in most of the data it is 
only by some 50% higher. Let us note that we arc using the old version of the UAl data 
[10] with rather arbitrary normalization, and we attempt to describe them only down to 
the lower limit of about 0.01 GeV^. In later publications of UAl [13] the BE ratio is 
shown to increase above 2 for lower Q^. Since, however, this increase is still a subject of 
controversy and anyway cannot be described by a Gaussian shape, we do not discuss it 
here. 

To explain why the BE ratio does not increase up to the value of 2, one may invoke 
some coherent component, but a more obvious effect (which also lowers the BE ratio) is 
the existence of longer living resonances. Pions coming from their decay are effectively 
"born" more than 10 fm from the collision point. Thus the Gaussian width parameter in 
a two-particle weight for these pions should by smaller by an order of magnitude, which 
allows practically to neglect their contribution to be BE effect in the experimentally 
accessible range. Therefore the Bialas - Krzywicki weights should be calculated taking 
into account only the permutations of momenta of pions produced directly, or resulting 
from the decay of the widest resonances. 

This is achieved easily if the procedure calculating weights is called before the decay 
of long-living resonances, i.e. in the same place, where the original LUBOEI procedure 
was called. We have rewritten correspondingly our program separating the procedure 
LWBOEI (called directly from JETSET, and calculating for each event a weight as a 
product of weights for positive, negative and neutral "direct" pions) from the master 
program (calculating distributions with and without weights). This procedure is available 
from authors as a FORTRAN file. The results are now (after rescaling the weights, as 
desribed above) quite similar to the data and may be brought to even better agreement 
by fitting the only free parameter - the Gaussian half-width of a two particle weight a. 
We discuss this comparison with data in more detail in the next section. 

3 Results and comparison with data 

We have generated 10^ events of pp minimum bias collisions at 630 GeV CM energy by 
the default version of the PYTHIA/ JETSET generator [6,11]. For each event the weight 
factor was calculated by taking the 4-momenta of "direct" pions of each sign, calculating 
for them a matrix of two - particle weights according to (2) with cr = 0.14 GeV, and 
then the weight w as a series (3) cut at X = 4 or 5. As already noted, the event weight 
is a product of weight factors for all three kinds of pions. 
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Since the charged pion multiphcity distribution (used to fit the default values of the 
parameters with all weights equal 1) is strongly affected by weights, we rescale the weight 
factors to restore the original average multiplicity. To this end we multiply the weights 
by a cF" factor, where n is the number of all "direct" pions of the event, and c and V 
are calculated from the comparison of the original and weighted multiplicity distribution. 
This is done by assuming that the original multiplicity distribution of " direct" pions may 
be well approximated by the negative binomial formula, i.e. that the NBD parameters n 
and 1/k are given by the experimental values of < n > and < n{n — l>/<n>^— 1. If 
with the weights we get a new average multiplicity < n' >, the original value may then 
be restored by rescaling the weights with 

^ ^ <n>{<n' > +k) 
<n' > {<n> +k) 

and 

^_ [1 + {1-V)<n'> /kY 
<w> 

where < w > is the average value of weights before rescahng. We have checked that this 
procedure restores indeed the original average multiplicity with accuracy of few percent. 
On the other hand, the BE ratios are little affected by rescaling (only the normalization, 
which is anyway mainly a matter of convention, changes by a few percent). 

In Fig.l we present the "BE ratio" (5) for pairs of positive pions as a function x — 
ln2{lGeV'^/Q'^) for the events from the original PYTHIA/JETSET generator (without 
weights) and from our prescription with series (3) cut aX K = A and K = 'b. In Fig.2 we 
show the "double ratio", i.e. the ratio of (5) for pairs of positive- and unlike sign pions 
for the same events. 




Fig.l. The "BE ratio" (5) for positive pious as a function of x = ln2{lGeV'^ /Q"^). 
Triangles, crosses and squares correspond to series (3) cut at K — (no BE effect), 
4 and 5, respectively. 
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Fig. 2. The "double ratio" of ratios (5) for ++ and +- pion pairs as a function of 
X = ln2{lGeV'^ / Q^) . Triangles, crosses and squares correspond to series (3) cut at 
K = (no BE effect), 4 and 5, respectively. 

We see that without weights both ratios are very close to one and depend weakly on 
(the dip in the double ratio at a; = 2.5, i.e. = O.lTGeV"^ is the reflection of in 
unlike sign pairs, and the wider dip at lower x comes from p). Our prescription produces 
a clear increase of both ratios at low Q^, and the difference between two choices of K^ax 
are almost neghgible. Thus we believe that cutting the series (3) at X = 5 we approxi- 
mate very well the results with full formula for the weights (1), which would require an 
unreasonably long computation time (even for supercomputers) when multiplicity exceeds 
20. 

In Fig. 3 and Fig.4 we compare our results obtained for < 5 and two values of a 
(0.14 and 0.1 GeV) with the UAl data [10] normalized as in Eq. (5). 
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Fig. 3. The "BE ratio" (5) for positive pions as a function of x = ln2{lGeV'^ / Q'^) . 
Crosses represent the UAl data [10], triangles and squares correspond to the series (3) 
cut at K — 5 with cr = 0.14 GeV and cr = 0.1 GeV, respectively. 
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Fig.4. The "double ratio" of ratios (5) for ++ and +- pion pairs as a function of 
X — ln2{lGeV^/Q'^). Crosses represent the UAl data [10], triangles and squares corre- 
spond to the series (3) cut at K — 5 with a — 0.14 GeV and a — 0.1 GeV, respectively. 

Let us stress here once more that we use on purpose an old version of the data, which 
seemed to be well described by a Gaussian shape of the ratio (5) for > O-OlGeV"^ 
{x < S) . Possible strong enhancement of this ratio for lower values of Q"^ discussed in 
more detail in the later UAl papers [13], which seems to signal a non-Gaussian shape, 
would require a modification of the shape of the two-particle weight Wij. Since here our 
purpose is merely to prove the reliability of our method we do not want to enter this 
problem. In most of the hadroproduction data such low values of are anyway not 
available and the data seem to be relatively well fitted by a Gaussian. 

We see that the data for the BE ratio up to x = 8 are bracketed by the results for two 
values of a, corresponding to the source radius of 1 and 1.4 fm. Wc do not attempt to fit 
a here more precisely, since our purpose is only to demonstrate the applicability of the 
method presented above. 

4 Conclusions and outlook 

We have shown that the weight method of implementing the Bose-Einstein interference 
effects in Monte Carlo generators may be applied effectively to describe the data. The 
prohibitive increase of computing time with multiplicity is avoided by an approximation, 
in which only the selected class of terms {K < K^ax in (3)) out of all n\ contributions 
is taken into account. We show that already for Kmax — 4 and K^ax — 5 the results are 
almost the same, which suggests that they approximate well those for the full series. 

The change of multiplicity distributions induced by weights is compensated by simple 
rescaling, equivalent to refitting of the original MC parameters. Using the simple one 
parameter Gaussian form of two-particle weight factor we describe reasonably well the 
UAl data. 
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There are many directions in which these results should be extended. One should 
check if the weight method can describe higher order BE effects and semi-inclusive ratios. 
The possibility of non-Gaussian and non-symmetric weight factors should be investigated. 
Other hadroproduction processes should be compared with pp collisions. We hope to learn 
soon whether our method is reliable enough to apply it confidently to the estimate of W 
mass shifts in the four-jet final states of the e+e" — W^W~ collisions. 
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